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ABSTRACT 


The topic of this thesis is the design of a toroidal prime mover. The primary 
design tool is a MATLAB program based on a standing wave analysis of thermoacoustics. 
A conventional prime mover that has been well studied, was used for validating the 
program. DeltaE, a program developed at Los Alamos National Laboratory for designing 
thermoacoustic engines, was also used as an auxiliary tool for validation test purposes. 
The validation suggested that the MATLAB program models the onset conditions of a 
conventional prime mover with fair accuracy. The program was then applied to a toroidal 
prime mover. A series of analyses were conducted on a toroidal prime mover to 
determine the influence of prime mover circumference, stack length, and duct radius. The 
final design calls for a 150 cm circumference, 2.54 cm duct radius, and 1.9 cm stack 
length. The contributions of this thesis are: (1) to develop and validate the MATLAB 
simulation program; (2) to provide a draft design for a toroidal prime mover for future 
research. 
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I. INTRODUCTION 


There have been many NPS theses involving several aspects of thermoacoustic heat 
transport. One area of particular concentration has been thermoacoustic prime movers 
[See Refs. 1-3, for example], A thermoacoustic prime mover is a class of thermoacoustic 
heat engine that converts thermal energy into useful work in the form of sound. The 
purpose of this thesis is to explore the design of a toroidal prime mover. 

A typical prime mover configuration is illustrated in Fig. 1. During operation, a 
temperature difference is imposed between the two heat exchangers, resulting in a 
temperature gradient in the stack. (The terms hot and cold are used to indicate that there 
is a temperature difference between these regions of the prime mover. They are not meant 
to imply anything about the absolute temperature involved). When the gradient is 
sufficiently high, acoustic oscillations spontaneously build up in the prime mover. At this 
point the prime mover is said to have reached onset. The conditions necessary for a prime 
mover to reach ortset are typically discussed in terms of which of two quantities is greater: 
the thermoacoustic gain of the stack, or the dissipation due to thermo-viscous wall losses 
and power delivery to an external load. The importance of each of these factors is 
determined by the geometry of the prime mover, e.g., where the stack is located relative to 
the nearest acoustic pressure antinode, etc. The positions of the components of a prime 
mover relative to the acoustic field are usually fixed by incorporating well defined acoustic 
boundary conditions into the design of the engine. For instance, placement of a rigid 
boundary in the engine insures that the acoustic velocity is zero at that point. These 
"built-in" boundary conditions serve as convenient starting points for computational 
analysis of the performance of an engine. 

A frequently asked question, when explaining how a prime mover generates sound 
(once the temperature gradient in the stack is sufficiently high) is "What is the origin of the 
low amplitude sound that gets amplified by the prime mover?" When answering, one 
usually resorts to postulation the existence of fluctuations or noise having a component 
with the right frequency and acoustic phase to be amplified by the stack. 
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Figure 1 - A typical prime mover configuration 
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In light of this discussion, a toroidal prime mover poses some interesting questions. 
As shown in Fig. 2, a toroidal prime mover has the cold end connected to the hot end via a 
duct that does not impose, at first glance, any obvious dominating boundary conditions 
that would force an antinode at any particular point. Indeed, if the stack and heat 
exchangers were not there and if the toroid had a uniform cross section, a standing wave 
should have no preferential orientation. In other words, there is no place that one knows 
a priori is going to be a pressure antinode, for example. However, one might argue that 
even with periodic boundary conditions, fluctuations still exist and there will be a 
component that matches the conditions necessary for amplification. 

Assuming that a toroidal prime mover can go into onset, a number of interesting 
questions arise. Will the standing wave (if it is an standing wave) adjust its spatial 
orientation to optimize the thermoacoustic gain? In other words, will the prime mover put 
the stack where it wants, where its performance is optimized? A rigid ended prime mover 
does not offer this possibility of self-optimization. Will deviations in measured stack 
locations from those predicted be indicative of unaccounted for processes such as acoustic 
streaming or turbulence? Another interesting question concerns multimode excitation. It 
has been previously shown that two modes can be excited in a prime mover 
simultaneously [Refs. 3 and 4], However, due to the presence of rigid boundaries, the 
stack placement is optimized for only one mode. A toroidal geometry with periodic 
boundary conditions offers the possibility that each mode could optimize itself. 

All these questions presuppose that the prime mover will go into onset. The 
presence of the stack will break the ambiguity of the standing wave's spatial orientation. 
The stack has the highest flow impedance in the toroid. Because this impedance will tend 
to favor low acoustic velocities, the perturbation due to the presence of the stack might 
force a velocity node at the stack. This condition would necessarily make any 
thermoacoustic gain zero, preventing the prime mover from ever going into onset. 
Further, the fundamental mode of a uniform cross section toroid (one without a stack), 
corresponds to the circumference being equal to one acoustic wavelength. There are two 
orthogonal degenerate modes that satisfy this condition. The presence of a nonuniformity 
will break this degeneracy. One mode will have a velocity node at the stack, the other will 
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have a pressure node at the stack. Neither of these modes will be supported by 
thermoacoustics, but linear combinations would be supported. Other nonuniformities in 
cross section, will also have an effect on the spatial distribution of the acoustic variables. 
Hence, by placing a constriction or expansion in the toroid, the previously mentioned 
velocity and pressure nodes will be displaced from the stack, allowing for the possibility of 
thermoacoustic gain. 

This discussion leads one to conclude that designing a functioning toroidal prime 
mover may be considerably more difficult than designing a conventional prime mover. 
The research described in this thesis addresses only the first part of the problem. That 
problem is to design a toroidal prime mover that, presupposing onset is possible, should 
reach onset at a temperature difference that is easily attainable in the laboratory. 

The approach taken is to use a simple formulation of the basic equations for 
thermoacoustics to predict the onset conditions for a toroidal prime mover. The program 
developed around this formulation will first be validated by comparing its predictions to 
data taken previously on a rigid ended prime mover [Ref. 1]. The predicted onset 
conditions will also be compared to those determined using DeltaE, a program developed 
at Los Alamos National Laboratory for designing thermoacoustic engines. After 
validation, the MATLAB program will be used to design a toroidal prime mover. 

Chapter II contains a brief review of the theory pertinent for this work. The 
validation of the program is discussed in Chapter III. Some experiments with a toroidal 
prime mover are presented in Chapter IV. An improved design is presented in Chapter V. 
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II. THEORY 


It has been shown previously that tracking the quality factor (Q) of a prime mover 
as a function of temperature difference across the stack is a convenient way to predict 
onset. This chapter contains a discussion of how to calculate the Q of a prime mover as a 
function of the applied temperature difference. The discussion is based on a standing 
wave analysis described in detail in Ref. 5 and draws heavily upon the derivation presented 
there. 

In addition, brief descriptions of the two different computer programs, one written 
in MATLAB [Ref. 6] by the author, and DeltaE [Ref. 7] will be given at the end of this 
chapter. 

A. STANDING WAVE ANALYSIS OF A THERMOACOUSTIC PRIME 
MOVER 

The derivation begins by making the connection between the power output of a 
prime mover and quality factor. Next, an expression is derived for the power output of a 
prime mover stack. Following this derivation, the thermal and viscous losses in the 
remainder of the prime mover are discussed. Finally, an expression is derived for the 
energy stored in the prime mover. 

1. Relationship Between Power Output And Quality Factor 

The reciprocal of the quality factor of a prime mover below onset can be expressed 
as 


l/Q = -W /g>Est 


( 1 ) 


The net power output W is the sum of the power output of the individual sections 
of the prime mover: 


W=W hot +W hothx +W pms +W coldhx +W cold . (2) 

The subscripts hot, hot hx, pms, cold hx, and cold refer to the hot end, the hot heat 
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exchanger, the prime mover stack, the cold heat exchanger and the cold end, respectively. 

Est is the acoustic energy stored in the prime mover, w is the angular frequency of the 
acoustic oscillation. At onset W- 0, so \IQ = 0. fV= 0 means that the power output of 

the stack exactly cancels the losses elsewhere in the prime mover, so the net power output 
is zero. 

2. Acoustic Power Output Of The Prime Mover Stack 

In Eq. (1) W pms is the power output of the stack. In general, W pms = \f x d W pms , 
where d W pms is the differential power output of a segment of the stack located between x 
and x + dx. It is found by taking the difference in the acoustic intensity at the two ends of 

the segment and multiplying the result by the cross-sectional area of the fluid in the stack 
■Spin* (the x direction is along the resonator axis, the direction of sound propagation). 

d oms = Spms[(pM^x+dx ~~ (pw)jc]. (3) 

The overbar indicates time average and the brackets () indicate average in the y direction. 
(The y direction is perpendicular to the resonator axis and the plate surface.) The stack 

plates are 21 thick and spaced by 2y 0 . Expanding (pu) x +dc and assuming that p is 
independent of>>, Eq. (3) can be written as 


dW pms = S pms ^-dx. (4) 

The acoustic variables can be expressed as p = p i e JWI and u = ii\e Jwt , where 

p\ =P^cos[<|)(x)], (5) 

and 


In these equations. Pa is defined to be the peak acoustic pressure amplitude at the rigid 
end, at x = 0. Changes in pressure amplitude due to changes in resonator cross section are 
neglected. Also, <|)(x) represents the phase of the standing wave at the location of the 
segment. If the propagation constant k is independent of position, <f>(x) = kx. However, in 

general. 


<K*) = fo KQdt, 


(7) 
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The term (1 + l/yo) arises from the requirement that volume velocity be conserved and the 
desire to relate the acoustic velocity in the stack to Pa ■ 

The time average of the product of two complex quantities such as pand (u) is 


p(u) = \Re\p x (u\)\ 


(~ denotes complex conjugate), Eq. (4) becomes 


dW pms =\S pm sRe(p l ^- + ^(u l ))dx 


The expression for dp\!dx and d(u\)/dx, are obtained from Swift [Ref. 8]: 


dp i N 

~Z~ = ~i — 7~\ u 0 
dx 1 -f v 


d(u i) _ j_ fk-f v VT m dp x j m 2 f Y-l .V 
dx (1 -a)(l + e*) T m dx 1+e/v^ 1 ’ 


The functions fk,f v and e s are given by 


tanh [(1 + j)yp/8ic] 

(i +j)y o/5* 

tanh[(l +y)> , o/5v] 

(1 +y>o/Sv 

p w c p 5ttanh[(l +y)yo/5^] 
p*Cj8,tanh[(l +j)l/5 s ] 


where the thermal and viscous penetration depths in the gas are defined as 


S& — J2K/p m Cp ®, 
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5 V = v /2(.i/p„,co. 


(16) 


The definition of the thermal penetration depth in the solid plate 5 S is similar to that of 5*. 
The Prandtl number a is given by 


a = c p \i/K. 


(17) 


Substituting Eqs. (5), (6), (10), and (11) into Eq. (9) gives the power output of the 


element 


d Wpms ~ ^SpmsRz 


f 


[ p m a l 

v2 


Jk-fv 1 + Hy p VT,„ . , (y- 1 ) ~ 

J\ ~ 7i r r rp tan cf)(jc) + 1 + —- 

\-f v (1 ~ CT )0 +Si) kT m YW \+z s J ) 


xcos 2 (x) ~j ~~^ sin 2 4)(x) 

* /v 




(18) 


The total acoustic power output of the stack is found by integrating Eq. (18) along the 
entire length of the stack. 

3. Thermal and Viscous Losses In The Remainder of The 
Prime Mover 

Thermal and viscous losses can be calculated from Eq. (18) with only slight 

modifications. The calculation is simplified by the absence of temperature gradients in the 
remainder of the prime mover, i.e. VT m = 0. 


dW hx = \S HX Re\ Pa 
2 I P nfl‘ 


1 + V cos2 ( t ) ^“ ~ sin 2< t ) ( x ) 


>dx. (19) 


Equation (19) also describes the losses in the cold end and hot end if / is set to zero and 
the appropriate form for S,f k , and/ v are used. For a circular tube of radius R, 

2J,[(/-l)/?/5,1 

' (J - 1 )(R/5 k )Jo [(/' - 1 )R/5 k ]' < 2 °) 

and 
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2Jx[(jj-\ym v \ 

Jv Q-i)(m v )j 0 w-i)R/d v ] 


(21) 


The thermal losses at the end caps are found by evaluating Eq. (19) at x=0 (the cold end ) 
and x=L (the hot end) and replacing .yo with L . 

4. Stored energy 

The energy stored in the prime mover is found by integrating the time-averaged 
acoustic energy density throughout the entire volume of the prime mover. The result is 



where x co id and x ho , refer to the length of the cold and hot end, respectively. 

B. MATLAB PROGRAM 

A MATLAB program was written to implement the standing wave analysis 
presented above to calculate 1 IQ vs. A T. In regions where there is no gradient 
(everywhere except in the stack), integrals are evaluated analytically. A linear gradient is 

assumed in the stack. The stack is divided into 10 isothermal segments. The 
temperature-dependent thermal-physical parameters and acoustic phase are evaluated for 
each segment and the work output computed. The total work output of the stack is the 
sum of the work output of each segment. In rigid-ended prime movers, the position of the 
elements of the engine relative to the standing wave are determined by boundary 
conditions. In the toroidal prime mover, the location of the elements relative to the 
standing wave is an adjustable parameter. 
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C. DELTAE 


DeltaE (Design Environment for Low Amplitude Thermoacoustic Engines) is a 
computer program written by Ward and Swift at Los Alamos National Laboratory [Ref.7], 
It solves Rott's wave equation for a geometry defined by the user. More information on 
DeltaE can be found on internet at http://rott.esa.lanl.gov. 
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m. VALIDATION 


The goal of this thesis is to design a toroidal thermoacoustic prime mover. As 
discussed in the previous chapter, the primary design tool is a MATLAB program based 
on the analysis presented in Ref. 5. The primary design goals are 1) to have an onset 
temperature difference that is easily reached experimentally and 2) to have dimensions that 
are practical from machining and experimental considerations. 

The purpose of this chapter is to explain the process used to validate the MATLAB 
program. The program was validated in two ways. 1) It was used to predict the 

performance of Lin's prime mover below onset, a system that has been well studied 
previously. 2) Predictions of onset temperature difference (A T 0 met) were compared to the 
predictions of DeltaE. 

A. LIN’S PRIME MOVER 

This prime mover, described in detail in Ref. 1 and Ref. 5, was used originally to 
investigate the low amplitude work output of a prime mover stack as a function of mean 
gas pressure and temperature difference across the prime mover stack. Figure 3 shows the 
prime mover configuration. 

Both the heat exchangers and prime mover stack have parallel-plate geometry. The 
work output of this prime mover stack is modeled by Eq. (18) in Chapter II. The 
thermo-viscous losses in heat exchangers are modeled using Eq. (19). The losses in the 
cold and hot ends are calculated using Eq. (19). Equation (20) and (21) were used in the 
calculation instead of Eq. (12) and Eq. (13). 

Figures 4 to 7 compare the MATLAB program results with the experimental results. 

The experimental data presented in these figures are composites of individual data sets 
recorded at different times. Figures 4 and 5 show how 1/Q changes with the AT for the 
prime mover filled with helium gas at a mean pressure of 376 kPa and 500 kPa. The 

computer simulation tends to underpredict the experimental onset values by about 20-30 
degrees. Figures 6 and 7 show how resonance frequency changes with AT. The 
difference between simulation and experiment is less than 6 Hz, or approximately 1%. 

One known limitation of the standing wave analysis is that it ignores boundary layer effect 
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Figure 3 - Lin’s prime mover configuration 
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0.03 


Lin's prime mover (pressure=376 kPa) 



Figure 4 - Graph of 1/Q vs. AT for Lin's prime mover filled 
with helium gas at a mean pressure of 376 kPa. 
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Lin's prime mover (pressure=376 kPa) 



Figure 6 - Graph of resonance frequency vs. A T for Lin’s prime mover 
filled with helium gas at a mean pressure of 376 kpa. 
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when calculating resonance frequency. However, the results are acceptable for the 
intended purpose. 

B. DELTAE PROGRAM SIMULATION 

The results of the DeltaE simulation of Lin's prime mover are shown in Figs. 8-10. 
The prime mover is filled with helium, at a mean pressure of 376 kpa. Figure 8 shows the 
heat input to hot heat exchanger as a function of the work output of the prime mover. 
Notice that there is a linear relationship. Extrapolation of a linear least squares fit to the 

points shows that the work output goes to zero at a heat input of approximately 44 W. 
Figure 9 shows the square of the acoustic pressure amplitude as a function of Q h . There 
is also a linear dependence between the two variables. Figure 10 shows AT as a function 
of Q h . AT has a minimum value of approximately 321 K, at a heat input of 
approximately 70 W. This minimum is taken to be the AT onse t. This value is indicated in 
Fig. 4 by the " * " symbol. The DeltaE result is in good agreement with both the 

measured result and the result of the MATLAB program. 

A similar procedure (see Figs. 11- 13) yields a ATonset of approximately 341 K for a 
mean gas pressure of 500 kPa. This result is indicated on Fig. 5 by the " * " symbol. This 

value is also in good agreement with the measured predicted results. Detailed input files 
for DeltaE are given in Appendix A and B. 
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QH 

Figure 9- Graph of the square of the acoustic pressure amplitude 
vs. heat input for Lin's prime mover filled with helium 
gas. The results are from DeltaE. The mean gas pressure 
is 376 kPa. 
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Figure 11 - Graph of heat input vs. work output for Lin's prime mover 
filled with helium gas. The results are from DeltaE. The 
mean gas pressure is 500 kPa. 
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Figure 12 - Graph of the square of the acoustic pressure amplitude 
vs. heat input for Lin's prime mover filled with helium gas 
The results are from DeltaE. The mean gas pressure is 
500 kPa 
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IV. TOROIDAL PRIME MOVER 


Having demonstrated the ability to predict the onset conditions for Lin's prime 
mover with the MATLAB program, the next step is to analyze a toroidal prime mover. 
One possible configuration is shown in Fig. 14. Heat is supplied to the hot heat exchanger 
and heat is removed from the cold heat exchanger, forcing a temperature gradient across 
the stack as usual. A rigid-ended prime mover would simply have rigidly-capped hot and 
cold (isothermal) ducts attached to the heat exchanger. However, in the toroidal 
geometry the hot and cold ducts are connected together, resulting in a temperature 
discontinuity at the junction. Of course, this discontinuity would be smoothed out by 
thermal conduction, so that there would be a gradual, continuous temperature change 
along the duct. This gradient will have some thermoacoustic effect that must be taken into 
account. More importantly, if the duct has a non-negligible thermal conductivity, a 
significant amount of heat will be conducted through it, perhaps requiring an excessive 
heat input to reach onset. 

To avoid the detrimental conduction through the duct, a toroidal prime mover might 
have two insulator sections, one attached to each of the hot and cold heat exchangers. 
The use of insulators will allow the duct to be held at ambient temperature while at the 
same time having little conduction through it. The disadvantage of the insulators is that 
the temperature gradient across them can be high. However, since the surface area is 
much smaller than that of the stack, any thermoacoustic effect will be small. 

Due to unfortunate time constraints, attemps to build a toroidal prime mover were 
limited to off-the-shelf components and those already on hand. This combination resulted 
in a compromised design. These compromises were understood from the outset, yet it 
was decided to complete the construction as a learning step toward a functioning toroidal 
prime mover. 

This chapter contains a description of the toroidal prime mover, analysis of the 
prime mover based on the MATLAB program, and a description of measurements 
performed, along with the results. 
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Figure 14 - Toroidal Prime Mover Configuration 
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A. TOROIDAL PRIME MOVER DESIGN 


Figure 14 shows an illustration of the toroidal prime mover. This section contains 
specific details. The hot and cold heat exchangers and the stack and stack holder are 

those used by Castro [Ref. 2], More complete details can be found there. 

1. Ambient Duct Assembly 

The ambient duct assembly is made from four sections of bent, galvanized steel 
conduit. Each section is 4.4 cm O.D., 4 cm I.D., and 48 cm long. The tubes are joined 

together with 2.35 cm long brass collars and soldered. 

2. Insulator 

The stainless steel insulators are sections of tube, 4 cm long, 3.85 cm I.D. and 4.0 
cm O.D. They are located between the ambient duct assembly and hot and cold heat 
exchangers to reduce heat conduction to the duct. The stainless steel tube is soldered to 
two stainless steel flanges. One flange is bolted to the ambient duct, the other to the heat 

exchanger and stack holder. 

3. Hot Heat Exchanger Assembly 

The assembly consists of hot heat exchanger, hot heat exchanger holder and hot heat 
exchanger heater collar. The copper heat exchanger is 3.924 cm in diameter and 0.257 cm 
thick. Each plate within the heat exchanger is 0.0292 cm thick. The plate gap is 0.0584 
cm. The heat exchangers are soldered into the heat exchanger holder. The copper heat 
exchanger holder is 7.366 cm O.D. and 0.953 cm thick. There is an O-ring groove in the 
side facing the stack. 

The hot heat exchanger heater collar is made of 0.953 cm thick copper and bolted 
around the hot heat exchanger holder. The collar is roughly square in cross section, 8.2 
cm long on edge. It has a 7.37 cm I.D. hole cut in it to fit over the heat exchanger holder. 
The collar is split and bolted together around the holder. There are two 0.625 cm 
diameter holes drilled in each piece to hold two heater cartridges, used to maintain a 

constant hot heat exchanger temperature. 

4.Stack Assembly 

The prime mover stack assembly consists of the stack and the stack holder. The 
stack is made of 0.01 cm thick polyester film. The film is rolled spirally around a phenolic 
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core 0.635 cm in diameter and 2.529 cm long. The polyester film is separated by 0.851 

mm diameter fishing line. The stack is 2.553 cm long, 3.772 cm O.D. and has a total 
surface area of 583 cm 2 . 

The 304 stainless steel stack holder is 2.576 cm long and 3.772 cm I D. (The stack 
holder is 0.22 mm longer than the stack creating a 0.11 mm gap on each side of the stack 
inside the holder). The stack holder has eight 8-32 machine screw holes in each flange to 
attach the insulator and hot heat exchanger on one side and the cold heat exchanger and 
insulator on the other. 

5.Cold Heat Exchanger Assembly 

The assembly consists of cold heat exchanger, cold heat exchanger holder and cold 
finger. In order to impose a temperature gradient across the stack, a 65.5 cm long, 0.7 cm 
thick, 10 cm wide copper plate (or cold finger) is bolted around the cold heat exchanger 
holder and immersed in liquid nitrogen. The cold heat exchanger assembly is similar to the 
hot heat exchanger assembly, except that the cold heat exchanger holder has no O-ring 
groove. The cold heat exchanger is sealed by using a 0.03 cm thick teflon gasket in the 
stack side of the cold heat exchanger holder and a 0.17 cm diameter lead wire diameter 
gasket on the insulator side. 

B. MATLAB PROGRAM SIMULATION 

In order to determine the best operating mean pressure, resulting in the lowest 
AT ' onS e t , the toroidal prime mover has modeled with the MATLAB program. As was 
mentioned earlier, designing a functioning toroidal prime mover is considerably more 

difficult than designing a conventional prime mover. For example, in a conventional prime 
mover, placement of a rigid boundary in the engine insures that the acoustic velocity is 
zero at that point. These boundary conditions serve as a convenient starting point for 
computational analysis of the performance of an engine. Unlike the conventional prime 
mover, the toroidal prime mover has no rigid boundary. Therefore, the location of the 

stack relative to the nearest pressure antinode is an adjustable parameter. 

Figure 15 shows a 3-D plot of A Tonset versus mean gas pressure of neon and x s , 
where x s is the distance from the location of the assumed pressure antinode to the hot end 
of the stack. The shape of the surface indicates the general dependence of A T onsel on 
those two variables. It is seen that the value of x s yielding the lowest A Tonset is 21 cm. 
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onset temperature difference(K) 



207 cm toroidal prime mover(pms length=2.529cm rad=1.905 cm) 



Figure 15 - Graph of A T omet vs. mean pressure and stack location. 
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Figure 16 shows AT onset versus mean gas pressure for x s =21 cm. It is seen that the 
lowest value of AT onset is approximately 190 K at a pressure of 100 kPa. Acknowledging 
the tendency of the program to under predict A T onsel , a reasonable value for A T omet is 
approximately 190-210 K. The largest possible A T attainable with liquid nitrogen and a 

hot end at nominal room temperature is approximately 220 K. Therefore, it is seen that 
the expected onset conditions are right at the limits of what is easily achieved. 

C. EXPERIMENTAL PROCEDURE 

The toroidal prime mover was pumped down using a vacumn pump to 
approximately 10 kPa and filled back with neon to approximately 100 kPa. This 
procedure is performed three times to ensure the gas inside the prime mover is 
substantially pure neon. Next, the dewar was filled with liquid nitrogen, cooling the cold 
heat exchanger. As the cold heat exchanger cooled down, heat was supplied to the hot 
heat exchanger via the heater cartridges to maintain the hot end temperature at 
approximately 293 K. The cold heat exchanger was surrounded by thermal insulation to 
reduce the heat leak to the surroundings. The temperature of the cold heat exchanger was 
measured with a type J thermocouple clamped between the cold heat exchanger holder 
and the cold finger collar. The coldest temperature achieved was approximately 110 K. 

Temperature control of the hot heat exchanger is achieved by a digital temperature 
controller (Omega CN 9000 temperature controller), heater cartridges, and a power 
supply amplifier (Techron 5530 power supply amplifier). The reference temperature for 
the temperature controller is taken from a type J thermocouple placed between the hot 
heat exchanger holder and the heater collar. The controller sends a 10 volt DC pulse to 
the mono input in the power supply amplifier. The duty cycle of the phase is proportional 
to the difference between the set point temperature and the thermocouple reading. Power 

to the heater catridges is provided by the output of the power supply amplifier. The hot 
end temperature was maintained to within approximately ±5 K of the set point 
temperature, initially 293 K. 

The pressure inside the toroidal prime mover is measured with a piezoresistive 
transducer with a sensitivity of 17.05 mv/psi. The transducer is mounted in the wall of the 
ambient duct, approximately 40 cm from the hot end of the stack. The transducer 
measures the mean pressure as a DC signal and the acoustic pressure amplitude as an AC 
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207 cm toroidal prime mover prime mover 



Figure 16 - Graph of h.T 0 „ se t vs. mean pressure of neon gas for x s = 21 cm 
as predicted with the MATLAB program. 

























signal, simultaneously. The transducer output is sent to a DC coupled preamplifier with a 
gain of 10. 

D. EXPERIMENTAL RESULTS 

The toroidal prime mover failed to reach onset at a mean pressure of 100 kPa and 
a AT of approximately 180 K (cold end temperature = 110 K, hot end temperature =293 
K). At this point the hot end temperature was increased to approximately 330 K, and the 

mean pressure varied over a range from approximately 20 kPa to 140 kPa. The prime 
mover still failed to go into onset. 

There are two possible explanations: 1) The applied AT was not large enough, or 2) 
The prime mover would never go into onset for reasons discussed earlier dealing with how 
the stack might influence the orientation of the standing wave. 

The null result suggests two possible paths for further study. 1) Design a prime 
mover with a lower A T onset This path is the subject of Chapter V. 2) Modify the existing 
toroidal prime mover to make it go into onset, to test the computer model. 

It was decided to remove approximately half the ambient duct from the toroidal 
prime mover and rigidly cap the remaining sections. This operation converted the toroid 
to a half-toroid. The new prime mover configuration is similar to a conventional prime 
mover except for the insulator sections and the two end ducts having the same 
temperature. 

The duct was cut and capped such that the hot end of the stack was 20 cm from the 

rigid end. The total length of the half-toroid was 108 cm. Figure 17 shows the output of 
the MATLAB program for the half-toroid configuration. The predicted AT onse , is 
approximately 117 K at a mean pressure of 100 kPa. 

The prime mover was filled with neon gas to a pressure of 100 kPa and the cold 

end was cooled down with liquid nitrogen, as explained previously. Onset was observed 
at a AT of approximately 110 K. This measured result agrees with the predicted value to 
within 7 K. The agreement gives some indication of how much the prediction for the 
whole-toroid might be in error. 
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temperature difference (K) 


half-toroid prime mover 


Pressure (kPa) 

Figure 17 - Graph of AT onset vs. mean pressure of neon gas for the 
half-toroid prime mover as predicted with the MATLAB 

program. The "o" represents the measured onset point. 
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V. IMPROVED DESIGN 


A. OPTIMIZING PARAMETERS 

The null result of the attempt to reach onset with the toroidal prime mover discussed 

in Chapter IV, leaves the basic question "Will a toroidal prime mover work?" unanswered. 
One obvious step is to redesign the prime mover to lower its predicted A Tonset to a value 
easily attainable. This redesign is the topic of this chapter. 

In considering the design, it was decided to assume the same stack and heat 
exchanger plate spacing as described in Chapter IV. The heat exchanger length is also the 

same. Three parameters were investigated: 1) the length of the stack, 2) the radius of the 

duct, and 3) the circumference of the toroid. It should be pointed out that A T omet for 
some of the parameters shown exceeds 216 K (the value of AT possible between nominal 

room temperature, 293 K, and liquid nitrogen temperature, 77 K). Therefore, for these 
plots we have set the hot heat exchanger temperature equal to 373 K. The temperature of 
the ambient duct is is 293 K. Thus, there is a 100 K temperature difference across the hot 

end insulator. This temperature difference is in the wrong direction, thermoacoustically. 
This will tend to make the predicted A Tonset higher than might be expected. Also, the 
thermophysical parameters of the gas will be somewhat different. However, these 

considerations do not significantly alter the conclusions draw from the results. 

Figure 18 shows the effect of changing the length of the stack. The graph shows 
A T 0 „ se t versus mean pressure of neon for x s = 21 cm. The symbols are for three different 
stack lengths: 3.175 cm, 2.529 cm, and 1.905 cm. The circumference is 207 cm. It is 

seen that the shorter stack gives the lowest value of A Tonset- There are two reasons for 
this result. First, for a fixed AT, a shorter stack results in a larger temperature gradient. 

Also, the shorter stack has less viscous loss. 

Figure 19 shows the effect of changing the radius of the prime mover. Again, the 

graph shows A T onset versus gas pressure at x s = 21cm and a 1.905 cm long stack. Three 
different radii are shown: 2.54 cm, 1.905 cm, and 1.27 cm. The lowest AT onset is occurs 

for the largest radius and is approximately 133 K. The reason for the decrease with 

increasing radius is that 1/Q is proportional to the ratio of the volume of gas contributing 
to losses to the volume of gas storing energy. This ratio is roughly proportional to 5/R. 
Finally, the effect of changing the circumference is investigated. As discussed in 

Chapter I, any nonuniformities in the cross section of the toroid may affect the spatial 
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orientation of the standing wave. Therefore, it is important to have as nearly a uniform 
cross section as possible. One way to make a uniform cross section is to machine a toroid 
out of a single piece of stock. A 207 cm circumference is perhaps too large to do this 
easily. Therefore, the effect of decreasing the circumference was investigated. In 
decreasing the circumference from 207 cm to 150 cm, the resonance frequency increases. 

Therefore, the position of the stack relative to the standing wave changes. In order to 

keep the stack at the right place, the value of x s must be decreased by a factor 150/207 to 
keep kx s the same in the two cases. Therefore, x s for the 150 cm circumference resonator 

is 15 cm. The results of this simple calculation are verified in Fig. 20 which is a 3-D plot 

for the 150 cm prime mover. It is seen that the minimum value of A T onset (152 K) occurs 
at x s approximately equal to 15 cm. Figure 21 shows A T onse i versus gas pressure for a 

1.905 cm long stack and 2.54 cm radius toroid. Two circumferences are shown: 207 cm 
and 150 cm. It is seen that decreasing the circumference increase A Tonset to approximately 
152 K. It also decreases the optimum mean pressure to approximately 70 kPa. 

Decreasing the circumference more increases AT onsel further. The 150 cm circumference 
is a reasonable compromise between the need for a low A T onset and having a prime mover 

size small enough for machining that is easily from a single piece of stock. 

B. SUMMARY 

The main finding of this thesis is a preliminary design of a toroidal prime mover. 

The results of a first-step design indicates that a useful geometry has a 1.9 cm long stack 
and 2.5 cm radius. AT on set for this configuration is predicted to be approximately 150 K at 
70 kPa mean pressure of neon. 
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The effect of changing the length of stack 
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Figure 18 - Graph of A T on set vs. mean gas pressure for the toroidal 
prime mover for three stack lengths. 









































The effect of changing the radius of stack 
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Figure 19 - Graph of A Tonset vs. mean gas pressure for the toroidal 
prime mover for three different duct radii. 
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Figure 21 - Graph of A T onsel vs. mean gas pressure for the toroidal 
prime movers for two different circumferences. 


























APPENDIX A. INPUT FILE FOR LIN'S PRIME MOVER 


TITLE Lin's Thermoacoustic Prime Mover Input File 
BEGIN Initial 0 

3.76E+05 a Mean P Pa 
591. b Freq. Hz G 
400 c T-beg K 
3.76E+04 d lpl@0 Pa G 
0.000 e Ph(p)0 deg 
0.000 fiUI@0 m A 3/s 
0.000 g Ph(U)0 deg 
helium Gas type 
copper Solid type 

ENDCAP Hot End 1 

1.146E-03 a Area m A 2 

sameas 0 Gas type 
copper Solid type 

ISODUCT Hot Duct 2 
sameas la a Area m A 2 
0.12 b Perim m 
0.055 c Length m 

sameas 0 Gas type 
copper Solid type 

HXFRST Hot HX 3 

sameas la a Area m A 2 
0.697 b GasA/A 
7.620E-03 c Length m 
5.200E-04 d yO m 
250.0 e Heatln W 
500.0 fEst-T K (t) 
sameas 0 Gas type 
nickel Solid type 

STKSLAB Prime Mover Stack 4 
sameas la a Area m A 2 
0.733 b GasA/A 
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3.500E-02 c Length m 
3.850E-04 d yO m 
1.400E -04 e Lpl ate m 
helium Gas type 
stainless Solid type 

HXLAST Ambient HX 
sameas la a Area m A 2 
0.697 b GasA/A 
2.040E-02 c Length m 
5.2E-04 d yO m 
0.000 e Heatln W (t) 
293. fEst-T K (t) 
helium Gas type 
copper Solid type 

ISODUCT Cold Duct 


sameas 

5a a Area m A 2 

0.12 

b Perim m 

0.8795 

c Length m 

helium 

Gas type 

ideal 

Solid type 

ENDCAP Cold End 

sameas 

5a a Area m A 2 

helium 

Gas type 

copper 

Solid type 

HARDEND 8 

0.000 

a R(l/Z) 

0.000 

bI(l/Z) 

helium 

Gas type 

copper 

Solid type 




APPENDIX B. INPUT FILE FOR LIN’S PRIME MOVER 


TITLE Lin's Thermoacoustic Prime Mover Input File 
BEGIN Initial 0 

5E+05 a Mean P Pa 
591. bFreq. Hz G 
400 c T-beg K 
5E+04 d lpl@0 Pa G 
0.000 e Ph(p)0 deg 
0.000 flUI@0 m A 3/s 
0.000 g Ph(U)0 deg 
helium Gas type 
copper Solid type 

ENDCAP Hot End 1 

1.146E-03 a Area m A 2 
sameas 0 Gas type 
copper Solid type 

ISODUCT Hot Duct 2 

sameas la a Area m A 2 
0.12 bPerim m 
0.055 c Length m 

sameas 0 Gas type 
copper Solid type 

HXFRST HotHX 3 

sameas la a Area m A 2 
0.697 b GasA/A 
7.620E-03 c Length m 
5.200E-04 d yO m 
250.0 e Heatln W 
500.0 fEst-T K (t) 
sameas 0 Gas type 
nickel Solid type 

STKSLAB Prime Mover Stack 4 
sameas la a Area m A 2 
0.733 b GasA/A 
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3.500E-02 c Length m 
3.850E-04 d yO m 
1.400E-04 e Lplate m 

helium Gas type 
stainless Solid type 

HXLAST Ambient HX 
sameas la a Area m A 2 
0.697 b GasA/A 
2.040E-02 c Length m 
5.2E-04d yO m 
0.000 e Heatln W (t) 
293. fEst-T K (t) 

helium Gas type 
copper Solid type 

ISODUCT Cold Duct 6 

sameas 5a a Area m A 2 
0.12 bPerim m 

0.8795 c Length m 
helium Gas type 
ideal Solid type 

ENDCAP Cold End 7 

sameas 5a a Area m A 2 
helium Gas type 
copper Solid type 

HARDEND 8 8 

0.000 a R(l/Z) 

0.000 b I(l/Z) 
helium Gas type 
copper Solid type 
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